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Scaling violations are found in the phase-ordering two- 
dimensional Heisenberg [0(3)] model, which has non-singular 
topological textures, under dissipative non-conserved dynam- 
ics. Three separate length-scales are found: Lt characterizes 
the scale of individual textures, Ln characterizes the sepa- 
ration between textures, and Lc characterizes the distance 
between oppositely charged textures. 
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Dynamical scaling has had great success in describing 
phase-ordering systems with purely dissipative local dy- 
namics ffl. The scaling hypothesis is that the evolution 
of the system at late times can be characterized by a 
single growing length-scale, L(t), which scales the two- 
point correlations of the system. As yet, however, there 
is no general approach to determine whether a particu- 
lar system scales. Conserved spherical systems can be 
shown to violate scaling 0] , while most other vector sys- 
tems that have been investigated seem to scale fairly well. 
Exceptions are provided by one- and two- dimensional 
systems with weakly interacting non-singular topological 
textures. Textures have both an internal scale charac- 
terizing the size of an individual texture, and an exter- 
nal scale characterizing the texture separation. If these 
scales evolve with different growth laws then scaling will 
be violated. 

The simplest system with topological textures is the 
one-dimensional (ID) XY [0(2)] model, in which the tex- 
tures are windings and anti-windings of the XY phase 
along the system. These, and all other, textures are non- 
singular, without the singular core structure of, e.g., vor- 
tices or domain walls. In ID, anomalous growth laws 
for both non-conserved and conserved dynamics || in- 
dicate scaling violations [|| that have recently been ob- 
served and explained Q. In this paper, we explore two- 
dimensional (2D) systems with textures, which we find 
also violate dynamical scaling. Unlike the ID systems, 
however, no extended scaling description is found to ap- 
ply to the two-point correlations. 

2D Heisenberg [0(3)] systems support topological tex- 
tures — also variously called skyrmions, instantons, or 
"baby-skyrmions" . An isolated texture can be pic- 
tured as a stereographic projection of an order-parameter 
sphere onto the plane of the system [Q. The intrinsic 
scale of the texture is then proportional to the radius 
of the projected sphere. These 2D textures have re- 
cently been of particular interest: e.g., in cosmological 
skyrmionic strings §|, in quantum antiferromagnets 0, 



and also in particle physics In all of these cases, 2D 
textures provide weakly-interacting, localized, but non- 
singular, excitations. 

Textures have an associated topological charge which 
is quantized. After an arbitrary sign choice, isolated tex- 
tures have total charge +1, while isolated anti-textures 
have total charge —1. Static solutions consisting solely 
of texture (or solely of anti-texture) configurations, in a 
(hard-spin) non-linear sigma model, were discussed by 
Belavin and Polyakov ||. The solutions are notable be- 
cause the energy of the system is independent of the over- 
all scale and also of the locations of the individual tex- 
tures, each one of which contributes unit charge and 87r 
energy (using Eq. (|^), below). These quasi-stable mini- 
mal energy solutions demonstrate the weak interactions 
of 2D textures. On the other hand, systems with both 
textures and anti- textures are never static |]]. However, 
the dynamics and evolution of generic textured systems 
is relatively unexplored fi|l3], particularly with dissipa- 
tive dynamics. Indeed, the "unwinding" , or annihila- 
tion, of textures with anti-textures appears not to have 
been addressed in two-dimensions. In this paper we con- 
sider mixed systems of textures and anti-textures devel- 
oped by instantaneously quenching a disordered state to 
T = 0. The phase-ordering dynamics are taken to be 
local, purely dissipative dynamics. For these systems, 
our numerical simulations find a rich pattern of scaling 
violations. 

The non-conserved "model A" dynamics are 



dt${x, t) = -6H/5$, 



where 



H = 



d 2 x 



{V<t>f + V a 
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The lack of thermal noise in (|l]) is appropriate for a 
quench to T = 0, where the equilibrium state has long- 
range order. We simulate the dynamics numerically with 
both soft (Vo < oo) and hard (Vb = oo, or equiva- 
lently with a \<j)\ = 1 constraint ) spins. We consider 
systems on square lattices of sizes between 128 x 128 
and 512 x 512, with periodic boundary conditions and 
independent randomly-oriented unit-magnitude spins as 
initial conditions. We use at least 20 independent runs 
for each system, and errors indicated in the figures are 
extracted from the variations between runs. We use a 
simple Euler update with a fixed time scale St = 0.01 
(except for Vq = 1/2, where we used St = 0.1), and con- 
sider times up to t = 10000. The late time regime covers 
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nearly three decades of time, starting with times t > 10. 
This regime is unchanged with a smaller timestep. Sim- 
ilar results are obtained with different Vo; although the 
early time behavior, and hence the asymptotic texture 
density, changes if Vo is small enough. In all cases, how- 
ever, the asymptotic growth laws remain the same. 

Previous numerical work has been carried out on these 
2D 0(3) systems, using Eqs. (0) and ((§), with both 
hard-spins (llj and with soft-spins (Vo — 10) |l2), al- 
though topological quantities were not measured. The 
hard-spin simulations, by Bray and Humayun, found an 
energy density consistent with e ~ t~ 2 / 3 , indicating a 
length-scale L ~ i 1 / 3 , and no dynamical scaling 111]]. The 
soft-spin simulations, by Toyoki, found L ~ ^o.42HT03 anc j 
dynamical scaling, however only the energy and spin-spin 
correlations at times t < 20 (in our units) were inves- 
tigated In this regime, we see strong early time 
transients. Our results are consistent with this previous 
work. However, we identify a distinct late time regime in 
which scaling is violated, and with growth law exponents 
that do not depend on Vo- 

The topological texture density for these systems is 
given by 

P (x) = [$-(djxd y $)]/4n, (3) 

which measures the local winding of the spin field around 
the unit sphere in order-parameter space. We show 
two snapshots of the texture density for part of a sin- 
gle 512 x 512 system in figures |l| and [2] For clarity, we 
only show texture density above the average magnitude 
in each figure. We identify the three evolving length- 
scales: Lt characterizes the scale of each texture, Ljv 
characterizes the separation between textures indepen- 
dent of charge, and Lc characterizes the separation be- 
tween oppositely charged textures. 

Dimensionally the equation of motion ([!]) determines a 
growing length-scale L ~ t 1 ! 2 , where we have suppressed 
the dimensioned kinetic coefficient. Growth laws differ- 
ent from this are only possible through the introduction 

— 1/2 

of other lengths: either the core scale £ ~ V , or 
the initial correlation length £o- We do not expect any 
£ dependence in this system since there are no singular 
defects. Any £o dependence, on the other hand, indicates 
a scaling violation, because then changing £o is not just 
equivalent to a shift in time So, for systems without 
singular defects, any growth law that differs from the di- 
mensionally naive one indicates a scaling violation. In 
addition, only positive powers of £o can enter into any 
growth law, yielding smaller growth exponents than 1/2, 
because shorter initial correlations should only decrease 
asymptotic length-scales. 

We first show the energy-density and the average mag- 
nitude of the texture density p in Fig. ||a. The en- 
ergy density e follows from Eq. (|^), and dimensionally 
scales as an inverse length-squared. The texture den- 
sity, from Eq. (H), has the same scaling. For spins of 
unit magnitude, (p) is a conserved quantity |jl3). How- 



ever is not, and it decreases along with e as tex- 
tures unwind with anti-textures. At late times e and 
fall off as £~°- 65±0 - 02 , indicating a growing length 
scale L N ~ f ^ ±om . [We fit times with t > 10 and 
extract the approximate errors from the variance of the 
exponents between different system sizes and different 
Vo-] This length scale, Ln, determines the overall texture 
density, and so characterizes the separation of individual 
textures. There is excellent quantitative agreement be- 
tween the different system-sizes, indicating that finite size 
effects only occur at later times, and qualitative agree- 
ment between the soft and hard-spin simulations at late 
times. For soft-potentials, the spins are not saturated 
at early times and this results in the initial increase of 
the average magnitude of the texture density seen for 
Vo = 1/2, and the shift in the asymptotic texture density 
at late times. 

We can also consider the portion of the energy con- 
tained in textures. Since an isolated hard-spin texture 
has an energy of 87r using (||) ||, we have plotted 87r(|p|) 
in Fig. ||a. The agreement with e at late times indicates 
that the asymptotic energy-density is wholly contained 
in the textures. The difference, e— (\p\) (not shown), de- 
cays very roughly like t~ l — and indicates a spin- wave 
like early time transient. This could explain why the 
strong early-time correction to scaling observed in e is 
much smaller in the texture density (|p|). 

While we find no way to scale two-point correlations 
of either spins or texture density, we can extract rele- 
vant length-scales from the evolution of the correlations. 
However, because of the lack of scaling, this analysis has 
significant systematic errors due to the evolving func- 
tional form of the correlations. In Fig. ^p, we show the 
position of the first zero of 

S(r,t) = <p(x)/>(x + r)). (4) 

This zero roughly characterizes the texture-antitexture 
separation, Lc- The growth saturates as Lc approaches 
the system size. We find Lc ~ ^o.4±o.i^ ^ over times 
t > 10 and excluding late times when finite size effects 
are apparent. 

There is a heuristic Energy-Scaling argument [Q for 
Lc ~ t 1 / 2 . We first identify the energy-density of 
the evolving system, e ~ 1/L 2 N . Independently, we 
calculate the rate of energy-density dissipation, dte = 
(dt<j) 8H/8cj>) = —{{d t (j>) 2 ), where we have used (|l|). Ap- 
plying the dynamics to ([2|), and neglecting the potential 
term since it is subdominant energetically [Q, we find 

dt4> — V 2 0. Inside a texture, V(j> ~ 1/Lt, while the 
Laplacian vanishes identically for pure texture solutions 
H . The natural length to enter in the second deriva- 
tive, then, is the texture-antitexture separation, Lc, so 
we expect V 2 ^ ~ X/LtLc- This will hold over the area 
of the texture (L 2 -,) for each texture (one per L 2 N ), result- 
ing in d t e ~ 1/ (L T Lc) 2 L\j L 2 N ~ \/L 2 c L 2 N . Comparing 
with e ~ 1/L 2 N , we find Lc ~ t 1 ' 2 . This growth law is 
consistent with the data. 
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What of the third length-scale, Lt? Since Lt is the 
texture scale, it characterizes the width of the peak cen- 
tered at r — of the two-point texture-density correla- 
tions. In Fig. ||c, we have plotted the width at half- 
maximum of 



AM)^<|p(xV(x + r)|)-<H) 2 



(5) 



Fitting the points for t > 10 we find L T ~ t°- 21±0 - 02 , with 
similar results from S(r, t). We can check this growth law 
by measuring (p 2 ). Since there is one texture, of area L T 
and with p 2 ~ 1/L T , in each area L 2 N , we expect (p 2 ) ~ 
1/(L 2 N L 2 T ). In Fig. §d, we measure (p 2 ) ~ r 1 - 02 * - 08 , 
by fitting t > 10. Using L 2 N ~ £°.65±o.02 f rom e, we 
have an independent estimate of Lt ~ i 019±0 02 . This is 
consistent, and we combine these estimates to find Lt ~ 

i0.20±0.02 

From the heuristic picture presented for Lc, small tex- 
tures have larger gradients and will annihilate earlier — 
leading to an increasing Lt- If we assume that the re- 
maining textures have not evolved, then their internal 
charge-density, 1/L T , is set by the initial fluctuations on 



the scale Ln, of order L^ /L 2 N . This implies Lt ~ L 



1/2 
J N 



Using our measurements for Lt and Ln, it is then likely 
that L N ~ t 1 / 3 and L T ~ i 1/6 (where (p 2 ) - t" 1 ). These 
exponents have been indicated by straight lines in Fig. 

I 

To summarize, we have found three characteristic 
length scales in the phase- ordering of the 2D non- 
conserved Heisenberg model. These length-scales exhibit 
different growth laws and so demonstrate the violation of 
dynamical scaling in this system. The growth laws are 
Lc ~ £°- 4 ± 01 describing the separation of textures and 
anti-textures, Ln ~ i°- 33 ± 01 characterizing the sepa- 
ration of textures, and Lt ~ ^o.20±o.02 describing the 
scale of individual textures. Soft and hard spin simula- 
tions give the same growth laws, despite individual tex- 
tures being destabilized by a finite Vb p[p^]. Similarly, 
the weak lattice instability |l3| does not seem significant. 
From our heuristic arguments, we believe that Lc ~ t 1 / 2 
exactly, and our data is consistent with Ln ~ t 1 / 3 ^ 3 



and Lt ~ t 1 / 6 ^ ■ The factors of the initial correlation 
length £o make up the dimensions of length needed for 
the scaling violations. A more convincing explanation for 
the growth laws awaits a detailed understanding of the 
texture-antitexture unwinding mechanism. 

To what extent does this picture of scaling violations 
carry over to other systems with non-singular topological 
textures? Certainly in one-dimension a similar picture 
holds H| . We believe that the observed scaling violations 
in ID and 2D are due to the weak interactions between 
textures. This being the case, we would expect similar 
scaling violations in conserved 2D 0(3) models, though 
these have not yet been explored to our knowledge. How- 
ever, because textures are strongly unstable for d > 2 
IP^I , it would be surprising if similar scaling violations 
were observed in higher dimensions. 
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FIG. 1. A snapshot of one-quarter of a 512 x 512 system 
with Vb = 10, at t = 13.1. Shown are the regions with texture 
density above the average magnitude, \p\ > (\p\). Black clus- 
ters have positive texture density, while outlined grey clusters 
have negative texture density. 



FIG. 2. A snapshot of the same system and region as in 
the last figure, at t — 27.9. Again, we show \p\ > (\p\). 
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FIG. 3. a) The energy density e (top plots) and the tex- 
ture density (\p\) (bottom plots) vs. time after the quench. 
The circles indicate indicate a system size Loo = 512 with 
Vo = 1/2, the squares indicate Loo — 512 with Vo = 10, and 
the asterisks indicate Loo = 128 with Vo = oo. The straight 
line has slope —2/3. The dashed line is 87r(|p|) from the 
Loo = 512 system with Vo = 10. b) The texture-antitexture 
separation Lc, characterized by the position of the first zero 
of S(r, t) vs. time after the quench. The symbols are the 
same, with stars indicating Loo = 256 with Vo = 256. The 
straight line has slope 1/2. c) The texture size Lt, character- 
ized by the width, at half-height, of the first peak of A(r,t) 
vs. time since the quench. The straight line has slope 1/6. d) 
The average square texture density vs time since the quench. 
The straight line has slope —1. 
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